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ABSTRACT 

We present new MMT spectroscopic observations of four dwarf galaxies represen- 
tative of a larger sample observed by the Sloan Digital Sky Survey (SDSS) and iden- 



tified by Peeples et al. (2008) as low-mass, high oxygen abundance outliers from the 



mass-metallicity relation. Peeples et al. (2008) showed that these four objects (with 



metallicity estimates of 8.5 < 12 -|- log(0/H) < 8.8) have oxygen abundance offsets of 
0.4-0.6 dex from the luminosity-metallicity relation. Our new observations extend 
the wavelength coverage to include the [O II] AA3726,3729 doublet, which adds lever- 
age in oxygen abundance estimates and allows measurements of N/O ratios. All four 
spectra are low excitation, with relatively high N/O ratios (N/O > 0.10), each of which 
tend to bias estimates based on strong emission lines toward high oxygen abundances. 
These spectra all fall in a regime where the "standard" strong line methods for metal- 
licity determinations are not well calibrated either empirically or by photoionization 
modeling. By comparing our spectra directly to photoionization models, we estimate 
oxygen abundances in the range of 7.9 < 12 -|- log(0/H) < 8.4, consistent with the 
scatter of the mass-metallicity relation. We discuss the physical nature of these galaxies 
that leads to their unusual spectra (and previous classification as outliers), finding their 
low excitation, elevated N/O, and strong Balmer absorption are consistent with the 
properties expected from galaxies evolving past the "Wolf-Rayet galaxy" phase. We 



compare our results to the "main" sample of Peeples et al. ( 2008 ) and conclude that 
they are outliers primarily due to enrichment of nitrogen relative to oxygen, and not 
due to unusually high oxygen abundances for their masses or luminosities. 

Subject headings: galaxies: abundances - galaxies: dwarf - galaxies: evolution 



1. INTRODUCTION 

There is a fundamental relationship between the mass of stars in a galaxy and its metallic- 
ity evolution (hereafter, the M-Z relation). Empirically, this has been observed as a luminosity- 
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metallicity relationship for low redshift dwarf galaxies (e.g., Lequeux et al. 1979, Skillman et al. 
1989 Lee et al.|2006b and references therein) and spiral galaxies (e.g., McCall et al.|[l985 Garnett 



&: Shields 1987; Zaritsky et al. 1994; Tremonti et al. 2004, and references therein). This robust 



relationship is observed over a range of 10 magnitudes in galaxy optical luminosity (e.g., Zaritsky 
et al. 1994 Tremonti et al. 2004; Lee et al. 2006b). In recent years, galaxies at higher redshifts 



have also shown a mass-metallicity or luminosity-metallicity relationship, and mounting evidence 



suggests this relationship evolves with time (e.g., Kobulnicky et al. 2003; Kobulnicky & Kewley 



2004 Shapley et al. 2004 Savaglio et al. 2005; Maiolino et al. 2008, and references therein), but 



see also Mannucci et al. (2010). Thus, this relationship provides both a very strong constraint on 



theories of galaxy evolution and a tool to better understand galaxies at higher redshifts. 



Tremonti et al. (2004, hereafter T04) derived luminosities, metallicities, and masses for ~53,000 



low redshift galaxies observed in the Sloan Digital Sky Survey (SDSS; |York et al.|2000 ) and convinc- 
ingly demonstrated that the basis of the empirically observed luminosity-metallicity relationship 
is an underlying association between stellar mass and metal abundance. The metallicities were 
estimated by fitting all observable strong emission lines and deriving a metallicity likelihood dis- 
tribution for each galaxy based on theoretical model fits calculated using a combination of stellar 
population synthesis and photoionization models. These models were then used with z-band lumi- 
nosities to form mass estimates. 

The physical driver for the M-Z relationship is still debated. Many studies favor supernova 
driven winds for the inability of a low-mass galaxy to retain its newly synthesized heavy elements. 



resulting in a lower effective yield with decreasing mass (e.g., Dekel & Silk 1986). Other observa- 



tional and theoretical studies are critical of this theory. For example, Dalcanton (2007) emphasizes 



the additional importance of star formation efficiency as outflows are an insufficient regulator in 
the absence of depressed star formation. Thus, a better understanding of the mass-metallicity 
relationship remains important. 

Although the M-Z relationship is well defined, it shows measurable scatter. Since observational 
error accounts for only half of the metallicity spread in the M-Z relation ( Cooper et al.|2008 >, one or 
more physical processes may be responsible for the remainder. Suggestions for the scatter include 
variations in the star formation history (e.g., recent starbursts, Contini et al. 2002), variations in 



stellar surface mass density ( Ellison et al.|2008 ), and variations in local galaxy density (e.g.. Cooper 



et al. 2008, and references therein). Motivated by the fact that outliers often provide key insights 
into the nature of physical relationships (e.g., Skillman et al. 1996), Peeples et al. (2008 hereafter 



P08) and Peeples et al. ( 2009 ) identified samples of outlier galaxies from the T04 study. 



P08 analyzed a sample of 41 high oxygen abundance, low-luminosity galaxy outliers from 
the M-Z relation of T04. Their "main" sample is comprised of 24 high abundance (8.95 < 12 + 
log(0/H) < 9.27), relatively low-mass (9.1 < log(M^/M0) < 9.9), low-luminosity (sub-L^; -19 < 
Mb < —17) dwarf galaxies. A redshift lower limit of z > 0.024 was imposed to ensure the inclusion 



- 3 - 



of the [O II] AA3726, 3729 emission doublet in the SDSS spectral coverag^ The data set was 
extended to include even lower mass galaxies, creating a second "very low mass" sample, with 7.4 
< log(M*/M0) < 9.0, -17 < Mb < -14, and 8.68 < 12 + log(0/H) < 9.12. In order to increase 
the sample to these lower masses, the redshift limit was dropped. All galaxies in the resulting "very 
low mass" sample have z < 0.024, and thus lack an [O II] AA3726, 3729 spectral measurement. The 
typical oxygen abundance offset in the 0/H - plane for a galaxy in the "main" sample is ~ 0.4 
dex, while the typical offset for the "very low mass" sample is ~ 0.6 dex (with offsets as large as ~ 
0.9 dex). 

P08 favored isolated or undisturbed systems with relatively low gas mass fractions nearing the 
end of their star formation activity as an explanation for the unexpectedly high oxygen abundances. 
Other possible causes for the high abundance outliers were considered and ruled out, such as 
discrepantly low luminosities for their masses or inaccurate metallicity calculations. 

Nebular oxygen abundances derived from the observations of strong emission lines and in the 
absence of direct measurements of the electron temperature are always subject to systematic effects 



(e.g., Kennicutt et al. 2003 and references therein). This is particularly true at the higher abun- 
dances typically found in spiral galaxies where several calibrations of methods based solely on strong 
lines result in systematically higher oxygen abundances when compared to oxygen abundances de- 
rived from direct measurements of the electron temperature (see, e.g., discussion in Bresolin]|2007 



and references therein). Importantly, Yin et al. (2007) and Perez- Montero &; Contini ( |2009 ) have 
pointed out that strong line calibrators that are based on the strength of the [N II] AA6548,6584 
emission lines are biased in the sense that large values of N/0 lead to overestimates of the oxygen 
abundance. At lower values of oxygen abundance, in general, the strong line methods show better 
agreement with the oxygen abundances derived from direct temperature measurements, but 



van 



Zee et al. ( |2006b ) have shown that there can be significant discrepancies at low values of excita- 
tion (A5007/A3727). Given the uncertainties in strong line oxygen abundance measurements, it is 
warranted to reinvestigate the conclusions of P08. 

The prospect of galaxy outliers from the M-Z relationship, and their consequences for galaxy 
evolution models, motivated the re-analysis of these objects. In this paper we discuss the previous 
SDSS and new MMT observations in § [2] and describe the analysis of the latter in § [3j In § [4] we 
use our new observations of the [O II] AA3726, 3729 emission lines to compare the properties of the 
outliers to the SDSS galaxies. Section[5]is dedicated to looking at several metallicity determinations 



and the appropriate applications, including the 03N2 method (§ 5.3.1), the N2 indicator (§ 5.3.2) 



and the R23 index (§|5.3.3). Our best estimates of the oxygen abundances, the nature of the objects 



in the "very low mass" sample that gives rise to their discrepant spectra, and an inspection of the 
nature of the objects in the P08 "main" sample are discussed in § [6| 



^ Seven galaxies with z < 0.024 passed the subsequent error and visual inspection cuts imposed by Peeples et al. 
( 2008 1 and so were kept in their "main" sample. 
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2. DATA 



2.1. Sample 

Four of the 17 metal-rich galaxies identified by P08 as "very low mass" (log(M^/M0) < 8.7) 
outliers from the M-Z relationship were selected for follow-up observations with the MMT (see 
Table [T]). These targets were chosen both for their significant departures from the mass-met allicity 
relationship and their availability during a single scheduled observing run. They have suggested high 
metallicities, as measured by T04, of 8.69 < 12+log(0/H) < 8.86. For SDSS J022628.28+010937.7, 
SDSS J024121.80+000329.2, SDSS J082639.19+253553.5, and SDSS J082633.77+252959.2 respec- 
tively (hereafter abbreviated as SDSS- plus the first six digits of the Right Ascension), these oxygen 
abundances lie 0.37, 0.40, 0.55, and 0.55 dex above (and well outside the ~ 0.1 dex scatter of) the 
luminosity- metallicity relationship histogram medians presented by P08. Like all of the objects in 
the P08 parent sample, these four galaxies are fairly isolatecQ low-redshift (0.0051 < z < 0.0227) 
dwarfs (—17 < < —15), with no obvious companions and somewhat depressed star formation 
rates (P08). 



2.2. SDSS Spectra 



The measurements made by T04 and used by P08 were derived from the publicly available 
SDS^ Data Release 4 ( |Adelman-McCarthy et al.| |2006[ ) data files referenced in Table [!} We 
used the SDSS pipeline reduced spectra (rather than performing our own 1-D extraction and 
reductions) to minimize differences between analyses. For a thorough description of the data 



^Note that J082639. 19+253553. 5 and J082633. 77+252959. 2 are members of the same group (Peeples, private 
communication, 2010). 

^http:/ /www. sdss.org/dr4/ 

Table 1. Present Sub-Sample of Low-Luminosity Outliers from Peeples et al. 2008 



Object 


SDSS Spectra 


RA 


DEC 


Ms 


log M* 


Redshift 


g-i- 


SDSS J022628.28+010937.7 


spSpec-51869-0406-561.fits 


36.6179 


1.16053 


-16.86 


7.92 


0.0051 


0.51 


SDSS J024121.80+000329.2 


spSpec-52177-0707-355.fits 


40.3408 


0.05813 


-16.56 


8.70 


0.0227 


0.46 


SDSS J082639.19+253553.5 


spSpec-52945-1586-164.fits 


126.6633 


25.59821 


-15.87 


8.59 


0.0078 


0.71 


SDSS J082633.77+252959.2 


spSpec-52945-1586-161.fits 


126.6407 


25.49979 


-15.23 


8.07 


0.0072 


0.55 



Note. — Publicly available SDSS spectra can be found at SDSS.org. The RA, DEC, Ms, log M*, and redshift values 
are taken from P08 (see their Table 1 for more details). The SDSS DR 7 provided g-r colors. 
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Fig. 1. Spectra of SDSS J022628. 28+010937.7. The two panels compare the optical H II region 
spectra from the SDSS and the MMT. Note that the SDSS spectral bandpass does not contain 
[Oil] A3727, but does encompass measurable [O II] AA7320,7330, as shown in the upper inset box. 
The bottom inset box is an enlarged view of the MMT spectrum to emphasize the dominance of 
[O II] over [O III]. This scale also reveals the significant [N II] AA6548,6584 strength and the lack 
of a [O III] A4363 detection. 
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Fig. 2.— Same as Figure[T} but the spectra are of SDSS J024121.80+000329.2. The SDSS spectrum 
is too noisy to detect the [O II] AA7320,7330 doublet, but displays the strong nitrogen lines. Again, 
note the strength of the [O II] A3727 emission relative to the [O III] emission in the MMT spectrum. 
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Fig. 3. — Same as Figure[T| but the spectra are of SDSS J082639. 19+253553. 5. The SDSS spectrum 
falls short of capturing the [O II] A3727 line and lacks the high signal to noise needed to detect the 
red [O II] lines. Below, in the MMT spectrum, the [O II] A3727 line strength exceeds that of the 
[O III] AA4959,5007 doublet. Also notice the strong [N II] AA6548,6584 lines. 
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Fig. 4.— Same as Figure[T} but the spectra are of SDSS J082633. 77+252959. 2. The MMT spectrum 
for this object has the most prominent [O II] A3727 hue relative to the [O III] doublet and, again, 
strong [N II] lines. Note that while [O III] A4363 appears to be measurable, its flux is less than the 
3 a uncertainty, and so was not used. 
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reduction refer to Adelman-McCarthy et al. (2006). While the median signal-to- noise (S/N) values 
for the SDSS spectra (12-47) meet the S/N (per pixel) > 8 requirement for reliable metallicity 



estimates (Kobulnicky et al. 1999), the wavelength coverage (3800-9200 A) excludes the O II 
A3727 emission line from their spectra. Although the red [O II] AA7320,7330 lines are included in 
the wavelength range, they are only detected in the spectrum of SDSS-022628. In Figures [T]|4j we 
have indicated the locations of the [O II] A3727 and [O II] AA7320,7330 emission lines. 

In §[4] and later in the paper, we compare galaxies from P08 to samples from the SDSS. In these 
comparisons, we use values of nebular emission line strengths, stellar absorption line strengths, star 
formation rates, and stellar masses based on SDSS spectra and photometry from the MPA-JHU 
data catalogu^ Stellar masses were determined based on fits to photometry following the work of 



Salim et al, 



( |2007[ f 



Kauffmann et al. (2003), and star formation rates (SFR) were based on Brinchmann et al. (2004) 

Objects with [O II], [O III], H/3, [N II], and Ha hue strengths less than 5 



and 

a were filtered out in order to increase the quality of this data set. 



2.3. MMT Spectra 

2.3.1. Observations 

New MMT observations were acquired in order to obtain improved S/N spectra and extended 
blue wavelength coverage including the [O II] A3727 line. The MMT data were taken with the Blue 
Channel spectrograph ( Schmidt et al.|[l989 ) on the UT date of 2008 November 1-2. Sky conditions 
were optimal with no cloud cover and sub-arcsecond seeing. A 500 line grating, 1" slit, and UV- 
36 blocking filter were used, yielding an approximate dispersion of 1.2 A per pixel, a full width 
at half maximum resolution of < 3 A, and a wavelength coverage of 3690-6790 A. Bias frames, 
flat-field lamp images, and sky flats were taken each night. The latter were primarily necessary 
due to significant differences between the chip illumination patterns of the sky and the MMT Top 
Box that houses the Blue Channel incandescent flat-field lamp. Multiple standard stars from |Qke 



(1990) with spectral energy distributions peaking in the blue and containing minimal absorption 
were observed throughout the night using a 5" slit over a range of airmasses. 

All four galaxies had strong central brightness peaks which were centered on the 1" x 180" 
slit. Three 900 second exposures (600 seconds for SDSS-086239) were made at a fixed position 
angle which approximated the parallactic angle at half the total integration time. This, in addition 
to observing the galaxies at an airmass below 1.5, served to minimize the wavelength-dependent 



light loss due to differential refraction (Filippenko 1982). A single slit position for each target 



was sufficient to characterize the global oxygen abundance, as metallicity gradients are small in 



*http:/ /www. mpa-garching.mpg.de/S DSS/ 

^MPA-JHU used a method similar to Sahm et al. (20071 to aperture correct their SFRs. 



low-mass galaxies (e.g., |Skillman et al.||1989[ [Kobulnicky fc Skillman||1996[ |l997[ |Lee et al.||2006a[ ). 
Finally, combined helium, argon, and neon arc lamps were observed at each pointing for accurate 
wavelength calibration. 



2.3.2. Data Reduction 



The MMT observations were processed using ISPEC2D (Moustakas & Kennicutt 2006), a 
long-slit spectroscopy data reduction package written in IDL. A master bias frame was created 
from > 20 zero second exposures by discarding the highest and lowest value at each pixel and 
taking the median. Master sky and dome flats were similarly constructed after normalizing the 
counts in the individual images. Those calibration files were then used to bias-subtract, flat-field, 
and illumination-correct the raw data frames. Dark current was measured to be an insignificant 
~ 1 e~ per pixel per hour and was not corrected for. 

Misalignment between the trace of the light in the dispersion direction and the orientation 
of the CCD detector was rectified via the mean trace of the standard stars for each night, pro- 
viding alignment to within a pixel across the detector. A two-dimensional sky subtraction was 
performed using individually selected sky apertures, followed by a wavelength calibration applied 
from the HeArNe comparison lamps taken at the same telescope pointing. Airmass dependent 
atmospheric extinction and reddening were corrected for using the standard Kitt Peak extinction 
curve dCrawford fc Barnes|p70l > . 



For each galaxy, the three sub-exposures were combined, eliminating cosmic rays in the process. 
The resulting images were then flux-calibrated using the sensitivity curve derived from the standard 
star observations taken throughout a given night. Finally, the trace fit to the strongest continuum 
source in the slit was used to extract the galaxy light within 15-30" apertures that encompassed 
> 99% of the spatially smooth and centrally peaked emission. Figures [T]|4] show the resulting one- 
dimensional spectra (with median S/N values > 30) in comparison to the SDSS spectra. Inset 
windows with a narrower spectral range emphasize the blue emission lines and the dominance of 
[O II] over [O III] in all four galaxies. 



3. ANALYSIS 

3.1. Emission Line Measurements 

Emission line strengths were measured for both the MMT and SDSS spectra using standard 
methods available within IRAFj^ In particular, the SPLOT routine was used to analyze the ex- 



®IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



- 11 - 



tracted one-dimensional spectra and fit Gaussian profiles to emission lines to determine their in- 
tegrated fluxes. The Ha and the adjacent [N II] lines were fitted simultaneously for an accurate 
deblending. 

Special attention was paid to the Balmer lines, which are located in troughs of significant 
underlying stellar absorption. Although the equivalent widths of the Ha emission lines were large 
enough that the underlying absorption was not a concern, this was not the case for H/3 and the 
bluer Balmer lines. We experimented with several different methods to correct the measurement 
of the H/3 emission flux for the effects of underlying absorption. At one extreme, maximizing the 
H/3 line intensity (by overestimating the underlying absorption) leads to a minimum Ha/H/3 value 
and a minimum estimate of the reddening. Because the underlying absorption is much broader 
than the emission and these features are both well resolved and relatively high signal-to-noise, the 
subtraction of the underlying absorption represents only a small contribution to the uncertainty in 
the reddening corrections (see Table [2]) . 

For our final analysis of the Balmer emission lines, we used the PA]>|^ analysis package to 
simultaneously fit the continuum, the Gaussian emission peak, and a broad, negative Lorentzian 
absorption feature. PAN uses a least-squares fit to minimize estimates the uncertainty 

using a "bootstrap" Monte Carlo error analysis, providing a reliable measurement of the fiux and 
associated uncertainty of the Balmer emission line. Figure[5]exhibits both the individual component 
and total fits determined by PAN, with minimal residuals validating the goodness of fit to our 
spectra. The emission line fiuxes are reported relative to H/3 in Table [2j and represent the multiple 
component fits for the H/3 and H7 emission lines and the single or deblended Gaussian profile fits 
from SPLOT for the rest of the emission lines. 

The errors of the flux measurements were approximated using 



where N is the number of pixels spanning the Gaussian profile fit to the narrow emission lines 
(typically 11). The rms noise in the continuum was taken to be the average of the rms on each 
side of an emission line. This error approximation is valid for weak lines whose uncertainty is 
dominated by error from the continuum subtraction. For cases where the fiux measurements were 
much stronger than the rms noise of the continuum, the error is dominated by flux calibration and 
de-reddening uncertainty. This is true for several lines in both the SDSS and MMT spectra of 
SDSS-022628, and for all the Ha lines (where a minimum uncertainty of 2% was assumed). Flux 
line strengths relative to H/3 and corresponding errors are listed in Table (2) As we were not able to 
detect [O III] A4363 and He II A4686 at the level of more than 3 a in any of the spectra, a fiux upper 
limit was estimated using Equation [T] We calculated upper limits on the electron temperatures 




(1) 



''pan was written by Rob Dimeo as a part of the Data Analysis and Visualization Environment, which is a 
software package developed at the NIST Center for Neutron Research for the reduction, visualization, and analysis 
of inelastic neutron scattering data, and was funded by the National Science Foundation. 
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Fig. 5. — region of the MMT spectra, as measured in PAN. The dotted gray hnes correspond to 
the Lorentzian absorption and Gaussian emission component fits. The red dashed hne displays the 
best overah fit to the original data (solid black line). The residual difference between the fit and 
data is shown in the bottom panel, with very little divergence near the emission peak, confirming 
a successful fit. The ensuing H/? emission fiux, accounting for Balmer absorption, can be found in 
Table [2 
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based on these A4363 fluxes, but none provided signiflcant constraints for abundance calculations. 
In some cases, the Ha/H/5 ratios for the SDSS and MMT spectra are significantly different. Using 
the MMT data for our reference ratio, the average percentage difference between the SDSS and 
MMT Ha/H/3 ratio is 15%, whereas the average [O III] A5007/II/3 percentage difference is 16%. 
Since the [O IIl]/H/3 ratio is less sensitive to flux calibrations, these differences between the two 
data sets are likely due to using long-slit (MMT spectra) versus circular fiber apertures (SDSS 
spectra). 



3.2. Reddening Corrections 



The wide range of observed wavelengths require fluxes to be corrected for extinction and 
reddening. Since the relative intensities of the Balmer lines are nearly independent of both density 
and temperature, they can be used to solve for the reddening. Assuming standard H II region 
characteristics (Te = 1.25 x 10^ K and Ue = 10^ cm'^), both the MMT and SDSS spectra were 



de-reddened using a Balmer decrement of 2.82 ( Hummer fc Storey] 1987) and the Cardelli et al. 
reddening law (with Ay = 3.1 E{B — V)). The extinction, ^i(A), from the Cardelli et al. 



1989 



1989) law was calculated using the York Extinction Solver (McCall 2004 
then derived the reddening value, E{B — V), using 



With these values we 



log 



I{Ha) 

Jim 



log 



F{Ha) 



+ 0.4 E{B - V) [Ai{Ha) - ^i(i^/3)]. 



(2) 



where F(Ha)/F(H/3) is the observed flux ratio and I(Ha)/I(H/3) is the de-reddened line intensity 



ratio using case B from Hummer & Storey (1987). Following Lee &: Skillman (2004) the reddening 



value can be converted to the logarithmic extinction at H[3 as 

c(i//3) = 1.43 E{B-V). 



(3) 



Original and de-reddened flux values for both the MMT and SDSS data are given in Table [2j 
where errors were propagated from those associated with the individual line measurements. All 
of the Ha/H/3 ratios are larger than 2.82, indicative of signiflcant extinction and reddening due 
to dust. Since the Galactic latitudes for the four galaxies are all large (ranging from 31° to 54°), 
little foreground extinction from Galactic dust is expected, and, indeed, the calculated values 
of E{B — V) (see Table [2]) are substantially greater than the foreground extinction determined by 
Schlegel et al. ( |1998 ). Our reddening corrections can be checked by comparing the corrected H7/H/3 
ratios with their theoretical values. As seen in Table [2| this ratio is consistent with the theoretical 
case B recombination ratio of 0.47 ( Hummer &: Storey||1987 ) in all four cases. This result implies 
accurate reddening corrections were made, which strengthens the argument for significant intrinsic 
extinction. 



http:/ /wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/YorkExtinctionSolver/ 
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COMPARISON OF SAMPLE GALAXIES TO THE SDSS 



Our new observations allow us to examine these Peeples et al. (2008) galaxies against larger 
samples drawn from the MPA-JHU SDSS database. Since the "very low mass" sample lacked [O II] 
measurements, we are particularly interested in how the strengths of these emission lines compare. 
Figure [6] illustrates properties of the four observed galaxies with respect to star-forming galaxies in 
the SDSS that fall in the "very low mass" sample range of 7 < log(M^/MQ) < 9. We have further 
restricted the SDSS comparison sample to those objects with s/n > 5 in the relevant lines. 

The top panel in Figure |6] shows a comparison of the excitation, measured by the [O II] /[O HI] 
ratio. The two more massive galaxies are clearly discrepantly strong in [O II] emission; these 
are very low excitation galaxies. With decreasing mass, the SDSS comparison sample becomes 
relatively sparse, but the two lower mass galaxies have comparatively low excitation. 

The second panel shows a comparison of the [N II] /[O II] ratio, primarily a function of the 
N/O abundance. Again, for the two more massive galaxies, there is a clear offset from the locus 
defined by the SDSS galaxies. While the [N II] /[O II] ratios are of comparable strength in the two 
lower mass galaxies, the trend in the SDSS sample is less clear. Extrapolating from the higher mass 
galaxies, the lower mass galaxies would be clearly discrepant, but they lie in the middle of a very 
sparse scatter in the diagram. Quantitatively, in a comparison with a least-squares fit to the SDSS 
compilation, the values of log(N/0) for our four observed galaxies are 0.2 — 0.4 dex higher than is 
typical for objects in the same stellar mass range. 

The lower two panels compare the equivalent widths of the underlying stellar absorption and 
the star formation rates derived from Ha. The four objects presented in this paper have H/3 
absorption equivalent widths on the high end of the SDSS distribution and star formation rates on 
the low side of the distribution. 

In sum, the four galaxies which we have observed are outliers in a number of properties. Given 
that oxygen abundances derived from strong emission lines are subject to a number of systematic 
uncertainties, it is clearly warranted to revisit their status as oxygen abundance outliers. 



OXYGEN ABUNDANCE DETERMINATIONS 



Relative to the SDSS spectra, our MMT spectra have the advantages of higher signal to 



noise and the inclusion of the blue [O II] A3727 line. Kniazev et al. (2003) showed from SDSS 



spectra that O^/H"*" ionic abundances can be determined reasonably well by observing the red 
[O II] AA7320,7330 lines (as a substitute for [O II] A3727); however, because they are auroral lines, 
their strong sensitivity to temperature can result in relatively high abundance uncertainties. In all 
four MMT spectra (Figures [iQ the preferred [O II] A3727 emission line strengths are noticeably 
stronger than the [O III] A5007 lines, underscoring the importance of accurately accounting for the 
contribution from the lower ionization state in determining oxygen abundances. 
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Fig. 6. — Properties of low-mass star- forming galaxies in the SDSS (gray dots) plotted against stellar 
mass. The four objects presented in this paper are shown as red dots in comparison. Interestingly, 
the [N II] /[O II] ratios are well above average for the mass range, and the [O II] /[O III] ratios 
show that these objects have significantly lower excitation than is typical. Additionally, their star 
formation rates are relatively low and their H/3 absorption equivalent widths are on the high end 
of the SDSS range. 
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Accurate "direct" oxygen abundance determinations from H II regions require a measurement 
of the electron temperature (typically via observation of the temperature sensitive [O III] A4363 
auroral line). However, as metallicity increases, cooling via metal lines becomes more efficient and 
the electron temperature decreases, making these intrinsically faint lines even more difficult to 
detect. Since we did not reliably detect [O III] A4363 in any of our targets, abundances must be 
estimated empirically or theoretically using relationships dependent upon relative strong line flux 
ratios. 



Alloin et al. 


(1979 


) and Pagel et al. 


(1979 



where oxygen abundance is related to one or more ratios of recombination and collisionally excited 
lines. Several other methods have since been developed and categorized as semi-empirical, empir- 
ical, or theoretical. Semi-empirical calibrations were determined using a combination of electron 
temperature measurements at low metallicity and photoionization models at high metallicity in 
correspondence with observational limitations. Empirical calibrations result from observations of 
H II regions with electron temperature measurements. However, the relatively small number of 
direct oxygen abundance determinations available are typically biased in the sense that they are 
based upon high-excitation H II regions only. In contrast, theoretical strong- line calibrations use ab 
initio photoionization models. One advantage to theoretical models is that they allow a wide range 
in ionization parameter in addition to input metallicity. However, these models rely on simplified 
assumptions of nebular properties and so do not yet provide entirely realistic representations of H II 
regions. Most troubling is that all theoretical strong line abundance determination methods over- 
predict abundances in the metal-rich (roughly solar metallicity and above) regime when compared 
to abundances determined from direct temperature measurements. 

Here we review the original oxygen abundances as determined by T04 using Bayesian models, 
the revised strong line oxygen abundances calculated by P08, strong line oxygen abundances from 
our new MMT spectra (with and without the addition of the blue [O II] emission), and the N/0 
relative abundance ratios. T04 found this type of comparison is valid for SDSS data by showing 
that analytic R23-metallicity relations roughly bracket the range of metallicities that they derive, 
concluding that their M-Z relationship is in line with previous strong-line calibrations. The cul- 
mination of these efforts is presented in Figure [7j which allows us to compare all of the oxygen 
abundance estimates with those expected from the M-Z relationship. Note that Figure [7] is shown 
for illustrative purposes only, and not to suggest which calibrator is more fundamentally correct. 



5.1. The T04 Oxygen Abundances 



P08 initially identified outliers from luminosities, gas-phase oxygen abundances, and stellar 
masses provided by T04. The T04 abundance calculations are based on a Bayesian statistical 



analysis of the strongest six emission lines using stellar population synthesis models from Bruzual 



fc Chariot] ( [2003^ and photoionization models from CLOUDY ( [Ferland et al.|1998[ ). Note that while 



the T04 abundances were used to identify the outliers (in a self-consistent manner), P08 derived 
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their own abundances, which are systematically lower than the T04 abundances. Yin et al. (2007 



hereafter Y07) compared the oxygen abundances from T04 with oxygen abundances derived from 
direct temperature measurements (using auroral lines) and found the T04 oxygen abundances to be 
systematically higher on average. Y07 further showed that the magnitude of this offset correlates 
well with the N/0 abundance ratio, and concluded that the offset is due to the assumption of a 



single N enrichment trend in the underlying Chariot & Longhetti (2001) modeling. 



Kewley Sz Ellison (2008 hereafter KE08) compared the M-Z relationships derived from ten 



different strong-line methods and found large systematic discrepancies between empirical and the- 
oretical calibrations. In the KE08 study, T04 metallicities tend to be in the middle of the range 
when low metallicities are considered (12 + log(0/H) ~ 8.5) and at the high end of the range 
when the high metallicities are considered (12 + log(0/H) ~ 9.0). KE08 provide relations for 
converting one metallicity scale to another, but the results of the Y07 analysis imply that a simple 
conversion is not always sufficient. That is, not only can the absolute abundances provided by the 
strong-line methods be systematically offset from the true nebular abundances, but, in some cases, 
the strong-line methods do not even accurately rank the abundances (see Table [s] and Figure [T] for 
different rankings among the 4 objects presented here). 



5.2. The P08 Oxygen Abundances 

P08 derived their own oxygen abundances for the outliers identified from the T04 sample. Since 
there was insufficient information to reproduce the results of the Bayesian statistical analysis of T04, 
P08 used two of the methods investigated by KE08. They used the [N 11] /[O II] ratio as calibrated 
by KE08 for the "main" sample (where [O II] A3727 was observed), and the [O III]/[N II] ratio (the 



"03N2" method) as calibrated by Pettini & Pagel (2004) for the remaining galaxies in the "very 



low mass" sample. P08 did not use the R23 diagnostic for their "main" sample because of concerns 
regarding its calibration at high metallicities (the presumed regime for their samples). P08 found 
that their oxygen abundances derived from the [N II] / [O II] ratio for the "main" sample agreed well 
with the T04 oxygen abundances. However, the 03N2 abundances were consistently lower than 
the T04 abundances by roughly 0.3 dex (in concordance with the offsets determined by KE08). 
The net effect is that the abundances derived by P08 are lower than those derived by T04. P08 
concluded that, regardless of the metallicity calibration used, their sample of 41 galaxies are true 
high-metallicity, low mass outliers from the mass-metallicity relation. Despite any offsets between 
the T04 and P08 data sets, both use a strong line calibration based partially on [N II] strength, 
which make abundance estimates liable to overestimates if the galaxies are nitrogen-enhanced (see 



5.3.2 and § 5.4). 
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5.3. Strong Line Methods Revisited 

We have chosen to re-compute the oxygen abundances for this sample using the MMT spectra 
with several strong line calibrations: the 03N2 method, the N2 method, the R23 index, and de- 
terminations using assumed temperatures. This serves to highlight differences amongst methods, 
outline the resulting range of possible abundances, and demonstrate that strong line calibrations 
are not appropriate for the present sample. These measurements are discussed below and compared 
in Figure [7} 



5.3.1. 'V3N2" Method 

The 03N2 method is one of two empirical calibrations used by P08 for the "main" sample, and 
the only calibration used for the "very low mass" sample since these objects lacked the necessary 



[O II] measurements. It was introduced by Pettini &; Pagel (2004, hereafter PP04) using empirical 



fits to strong line ratios from H II regions with "direct" oxygen abundances. Derived for the purpose 
of measuring metallicities in galaxies at high redshift, lines close in wavelength are used to mitigate 
the need for flux calibrations and reddening corrections. Direct metallicities were compared to 
the ratio of ([O III] A5007/H/3)/([N II] A6584/Ha) for a sample of 137 H II regions. See the 
resuhing 03N2 relationship of Equation 3 in PP04, where 03N2 = log(([0 III] A5007/H/3)/([N II] 
A6584/Ha)). 

P08 favored the 03N2 method with its high sensitivity to oxygen abundance. Indeed, 03N2 is 
a good calibration in the high-metallicity regime, where [N II] tends to saturate, but the strength of 
[O III] continues to decrease with increasing metallicity. Note that the choice by PP04 to use only 
objects with direct metallicities, in an effort to provide a more secure calibration, introduces a bias 
because low excitation spectra (i.e., relatively low values of A5007/A3727) are excluded from the 
calibration. Thus, the unintended consequence of the choice to limit the sample to objects which 
were perceived to have more accurate abundances has resulted in a biased sample. This biased 



sample gives the impression of a smaller scatter in the relationship than occurs in nature (cf., Yin 



et al. 2007). 



To confirm consistency with P08, we duplicated the 03N2 measurements for the SDSS spectra 
and calculated 03N2 from our MMT spectra. The results are listed in Table [3] and closely match 
the measurements by P08 using the same 03N2 methodology; they are lower than the values found 
by T04 by 0.2 - 0.3 dex. 

Y07 showed a large scatter in the comparison of abundances derived from the 03N2 calibration 
with oxygen abundances derived from the direct method. They suggest that the scatter could 
be due to excluding the ionization parameter in the 03N2 calibration. The importance of the 
ionization parameter strengthens in the low-metallicity regime where 03N2 is much less dependent 
on metallicity. Additionally, at low Z, N is thought to be a primary element, implying that the 
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N/0 ratio is independent of 0/H. 



Since the calibration of the 03N2 method does not account for possible variations in N/0 



for a given 0/H, inherent biases are possible. Perez-Montero &; Contini (2009, hereafter PMC09) 



found that a strong correlation exists between metallicity derived from the 03N2 parameter and 
the N/0 ratio, such that when the N/0 is enhanced, the 03N2 strong line calibration tends to 
overestimate the 0/H abundances. PMC09 corrected for this dependence by using 12 + log(0/H) 
versus 03N2 residuals to modify the 03N2 calibration (see Equation 8 in PMC09). Using the N/0 
calculations from § 5.4 and given in Table [41 we calculated the revised 03N2 abundances and list 



them in Table [sj This calibration lowers the 03N2 abundances only slightly (by an average of ~ 0.1 
dex) relative to PP04 03N2 estimates. In § 5.5 we look at the higher abundances from the 03N2 



method relative to relationships predicted by direct and photoionization model abundances. 



5.3.2. "N2" Indicator 



Denicolo et al. (2002, hereafter D02) proposed the use of the [N II] A6584/Ha ratio as a 
sensitive metallicity indicator. This is a promising estimator: like 03N2 it eliminates uncertainties 
due to reddening and flux calibrations, and functions well for metallicity values below the [N II] 
saturation level. The original N2 calibration developed by D02 (see Equation 2 in D02), where 
N2 = log([N II] A6584/Ha), was based on ~155 H II regions and used a least squares fit with an 
estimated uncertainty of ~0.2 dex. For our objects, this method produces oxygen abundances (8.61 
< 12 + log(0/H) < 8.70) in agreement with the P08 values. 

PP04 re-calibrated this parameter, primarily with electron temperature-based metallicities, 
characterizing the fit both linearly and with a third-order polynomial. The polynomial fit is given 
in Equation 2 of PP04, yielding oxygen abundances of 8.45 < 12 -|- log(0/H) < 8.56. The results 
for both calibrations are listed in Table |3l 

Note, again, as with the PP04 calibration of the 03N2 method, the exclusive use of calibrators 
determined with measured electron temperatures underestimates the scatter in the relationship, and 



biases against low-ionization H II regions. Yin et al. (2007) found similar results in their analysis 



of the N2 method as found for the 03N2 method. The uncertain roles of the ionization parameter 
and the relationship between N/0 and 0/H as a function of metallicity result in appreciable 
uncertainty and scatter in the N2 method. Expanding on this idea, PMC09 found an expectedly 
strong dependence of metallicities predicted by N2 on the N/0 ratio. PMC09 derived a calibration 



correcting for this effect (see their Equation 13). Using N/0 ratios from § 5.4 and tabulated 
in Table |4| we present the results of this calibration in Table [3j This method predicts oxygen 
abundances which are ~ 0.1 and ~ 0.2 dex smaller on average than those predicted by the PP04 
and D02 N2 methods respectively. Shifts to smaller oxygen abundances of this size can account for 
roughly half of the offset from the L-Z and M-Z relationships. 
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5.3.3. R23 Index 



Pagel et al.| ( |1979D promoted the use of the R23 index, R23 = ([O II] A3727 + [O III] AA4959, 
5007)/(H/?), as a good estimate of oxygen abundance in the absence of an electron temperature 
measurement. Because the optical [O II] and [O III] emission lines decrease at both high abundances 
(due to an increasing role of fine structure line cooling) and low abundances (due to the decrease 
in the relative number of oxygen atoms), the relationship is bi- valued and, therefore, potentially 
ambiguous. The turn-around in the relationship occurs at oxygen abundances near 12 + log(0/H) 
~ 8.4 where R23 reaches a maximum of ~10. However, the degeneracy between the two branches 
of solutions can, in most situations, be broken with an additional determinant. Also, since R23 is 
based on emission lines with significant separation in wavelength, accurate reddening corrections 
and uncertainties are important. Below, we discuss two R23 calibrations, one constructed from 
photoionization models, and another empirically based. Since P08 believed the sample to be high- 
metallicity, they did not consider the R23 diagnostic due to the potential of [O II] + [O III] to 
saturate and the lack of [O II] measurements in their "very low mass" sample. 



McGaugh (1991) created a calibration based on theoretical photoionization models using R23 
and the additional O32 index: O32 = ([O III] AA4959, 5007)/([O II] A3727). McGaugh developed 
this model using the photoionization code CLOUDY ( Ferland et al.|[r998 ) and zero-age H II region 
models, accounting for photoionization parameter variations. To discriminate between the two 



branches, van Zee et al. (1998), followed by others, advised using the ratio of ([N II] A6584)/([0 II] 





van Zee et al. ( 


1998 


iMcGaugh 


( 


1994 


) sugj 



and > 0.1 for high abundances, giving a rough distinction between lower and upper branches. Note, 



however, that McGaugh ( 1994) also points out (his Figure 3) that there is a dependence on the N/0 
ratio on the behavior of this discriminant. For the four galaxies in our sample, [N II] /[O II] ranges 
from 0.15 to 0.22 with an average uncertainty of ^^0.02, naively suggesting that they are upper 
branch objects near the turn-around region. We determined oxygen abundances using the analytic 



equations of the semi-empirical calibration from Kobulnicky et al. (1999), which have an estimated 



accuracy of ~0.15 dex ( [Kobulnicky k Kewley|[2004l > . The results for both upper and lower branch 
calculations are listed in Table [3} The upper branch values define an upper limit to the metallicity 
estimates and are roughly equal to the P08 values, while the lower branch values are approximately 
0.3 to 0.5 dex lower and would be consistent with the scatter in the luminosity-metallicity relation. 



Pilyugin et al. (2001a) and Pilyugin et al. (2001b) empirically calibrated the R23 index with 



electron temperature based H II region oxygen abundances. This relation was later refined by Pilyu- 



gin Sz Thuan (2005, hereafter PT05), where the resulting fit incorporates the excitation parameter 
(which accounts for the effect of the ionization parameter) P = ([O III] AA4959,5007/H/3)/R23. 
The upper and lower branches correspond to electron temperature based metallicities with 12 + 
log(0/H) > 8.25 and 12+ log(0/H) < 8.0 respectively (see Equations 22 and 24 in PT05). The 
PT05 model upper and lower branch values, both of which are several tenths of dex below the P08 
values and have an estimated accuracy of ~0.1 dex, are given in Table [3| However, note that all 
four of our objects lie outside of the calibrated region for the PT05 method (see PT05, Figure 12). 
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The [O III]/[0 II] ratio is sensitive to the ionization parameter, and important for characterizing 
the physical conditions in the H II region. Since the age of an H II region is hnked to the evolution 
of the ionization parameter and shape of the ionizing spectrum, photoionization models which 



assume a zero-age main sequence can lead to systematic errors, van Zee & Haynes (2006a) and 



van 



Zee et al. (2006b) discuss these effects and suggest they may be significant for H II regions with 
log(032) < —0.4, such as for three of the four objects in our sample (see Table[3]and §|6|. The small 
A5007/A3727 ratios in these spectra are indicative of a low ionization parameter, and correspond 
to very small P values (see Table [s]) that are inconsistent with our high values of log(R23) for the 
PT05 upper branch calibration. Similarly, our low P values fall outside the range of calibration for 
their lower branch. If we naively extrapolate their calibrations, these objects fall close to or within 
the turn-around region which has an abundance range of 8.0 < 12 -|- log(0/H) < 8.4. 

[N II] /Ha is often used as a second indicator of branch division. However, as it is less sensitive to 
metallicities and more responsive to ionization than [N II] /[O II], the division is less clearly defined. 
The latter ratio loses the ability to clearly break the degeneracy near the turn-around point of 
[N Il]/[0 II] ~ 0.1 (or 12 -|- log(0/H) ~ 8.4). With its dependence on both oxygen abundance 
and excitation, R23 seems like a superior metallicity indicator, but the ambiguity between branches 
requiring a secondary measurement like [N II] /Ha or [N II] /[O II] clearly render it less than ideal. 
The enhanced [N II] /[O II] values for the four objects investigated in this paper result in highly 
uncertain R23 metallicity estimates. Note, however, that the turn-around between the two branches 
happens at relatively low oxygen abundances, for all values of excitation, already indicating that 
the higher oxygen abundances reported by T04 and P08 are less likely. 

Finally, we note the possible utility of the red [O II] emission lines in low redshift SDSS spectra. 



Kniazev et al. (2003) showed that SDSS abundances estimated using the red [O II] AA7320,7330 
lines as a substitute for [O II] A3727 yield comparable results for spectra with direct electron 
temperature measurements. Since the red [O II] emission lines were detected in the SDSS-022628 
spectrum (Figure [T]) we were able to derive alternative R23 abundance estimates for that galaxy. By 
using the relative emissivities to estimate the relative strengths of the blue [O II] emission lines from 
the red [O II] emission lines, we found that the A5007/A3727 ratio is less than unity (suggesting 
relatively low ionization) over the range 1 x 10'^ K < Tg < 1.25 x 10^ K (for Ug = 100 cm~^). The 
corresponding R23 calibrated 0/H values are consistent with our MMT spectra findings, but, of 
course, the branch choice ambiguity remains. 



5.4. N/O Relative Abundances 

Since many of the methods of abundance determination discussed here rely on N abundance 
values either to calibrate 0/H values or to discriminate between bi- valued solutions, a confident 
measurement of the N abundance is important. Due to the relative insensitivity of the derived N/O 
ratio to electron temperature, it is possible to get a reliable estimate of this ratio in the absence 
of an electron temperature determination. We calculated the N/O values for our four galaxies (see 



Table 3. Metallicity Estimates for H II Regions Using Empirical Methods 



MMT Data 



Quantity 


SDSS-022628 


SDSS-024121 


SDSS-082639 


SDSS-082633 


03N2 


0.77±0.05 


0.45±0.08 


0.28±0.11 


0.66±0.13 


N2 


-0.71±0.04 


-0.63±0.06 


-0.58±0.06 


-0.69±0.09 


log R23 


0.60±0.02 


0.71±0.04 


0.64±0.06 


0.71±0.05 


log O32 


-0.20±0.04 


-0.64±0.08 


-0.74±0.12 


-0.50±0.11 


[NlI]/[OlI] 


0.22±0.01 


0.16±0.01 


0.20±0.02 


0.15±0.01 


[N II] /Ha 


0.20±0.01 


0.24±0.01 


0.26±0.02 


0.20±0.02 


P 


0.39±0.01 


0.19±0.01 


0.15±0.02 


0.24±0.03 


N/O 


0.16±.01 


0.12±0.01 


0.15±0.02 


O.lliO.Ol 


Pettini and Pagel (03N2) 


8.48 


8.59 


8.64 


8.52 


Perez-Montero and Contini (03N2) 


8.37 


8.52 


8.54 


8.47 


T\ 11'' 1 T~t ^ /TVTn\ 

Pettmi and Pagel (N2) 


8.45 


8.51 


8.56 


8.46 


Denicolo (N2) 


8.60 


8.66 


8.70 


8.61 


Perez-Montero and Contini (N2) 


8.30 


8.44 


8.42 


8.41 


Kewley and Dopita (N II/O II) 


8.83 


8.74 


8.80 


8.71 


McGaugh Upper (R23) 


8.77 


8.63 


8.70 


8.64 


McGaugh Lower (R23) 


7.81 


8.13 


8.06 


8.08 


Pilyugin Upper 


8.40 


8.03 


8.07 


8.11 


Pilyugin Low 


7.58 


7.43 


7.22 


7.55 


van Zee expected from M_b^ 


8.22 


8.17 


8.07 


7.97 






SDSS Data 






SDSS-022628 


SDSS-024121 


SDSS-082639 


SDSS-082633 


Pettini and Pagel (03N2) 


8.53 


8.61 


8.68 


8.50 


Pettini and Pagel (N2) 


8.48 


8.54 


8.59 


8.43 


Denicolo (N2) 


8.63 


8.69 


8.73 


8.58 


Tremonti (T04) 


8.82 


8.90 


8.86 


8.69 


Peeples 08 


8.51 


8.61 


8.66 


8.49 



Note. — Values for various methods of determining metallicity estimates from the strong lines of the 
MMT spectra are listed, with the 03N2 and N2 estimates for the SDSS spectra given below for comparison. 
Interestingly, the N2 method produces similar values for all 4 galaxies, and the 03N2 method reproduces 
similar values to those found by P08. 



^Metallicity estimates calculated using Equation [4] of this paper, as determined by 



van Zee et al. 



(2006b I. 
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(a) Luminosity-Metallicity Relation 



(b) Mass-Metallicity Relation 



Fig. 7. — The figure (a) shows 12 + log(0/H) vs. for our sample. The soHd hne is the Unear fit to dwarf 



galaxies found by Lee et al. ( 2006b ) for comparison, with the dashed lines representing the 1 a uncertainty. 



This relationship is similar to that of van Zee et al. ( 2006b ) , but provides an average dispersion for both the 



L-Z and M-Z relationships. The points plotted display calculated oxygen abundances for our sample. Each 
method is denoted with a different symbol, where objects are differentiated by color: green for SDSS-022628, 
purple for SDSS-024121, blue for SDSS-082639, and pink for SDSS-082633. The 'x' symbols denote the 
oxygen abundances calculated from emissivities assuming various electron temperatures for SDSS-082633. 
Figure (b) displays 12 + log(0/H) vs. log(M^,/M0). Again, the solid line depicts the linear relationship and 



1 a uncertainty of Lee et al. ( 2006b ) for mass. 
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Table [4]), assuming that N"*"/©"*" is approximately equivalent to that of N/0 (based on their similar 
ratio of ionization potentials, Vila Costas &: Edmunds||1993 ) and an electron temperature of 12,500 
± 2,500 K. 



The N/0 versus 0/H trend is well studied in galaxies of varying types. Vila Costas & Edmunds 



( 1993 ) presented a thorough overview of theoretical expectations and observations available at the 



time. A salient point is that N can be produced as both a primary and a secondary element and 
that the secondary component is expected to be delayed relative to oxygen and to dominate at 



high abundances. Garnett (1990) convincingly demonstrated that there can be significant scatter 



in N/0 at a given 0/H, and Izotov & Thuan (1999) confirmed that this scatter is significant for 
oxygen abundances above 12 + log(0/H) = 7.7. Garnett (1990) proposed that much of the scatter 



could be explained by the time delay between producing oxygen and secondary nitrogen. However, 



Henry et al. (2006) conclude that the scatter in N/0 could be produced by a variety of causes. 
As already indicated in § [4j the four objects in our sample have unusually high values of 



N/0. PMC09 and, later, Amorm et al. (2010) (in their study of the abundances of the "green 
pea" galaxies) present a modern view of N/0 versus 0/H assembled from a sample of 475 H II 
objects and the SDSS DR7 catalog of observations respectively. A clear trend is obvious, but with 
scatter in log(N/0) on the order of 0.5 dex or more for a given 0/H. An inspection of Figure 2 
from Amorm et al. (2010) shows that our four galaxies, which range between —0.8 < log(N/0) 
< —1.0, correspond to values of 12 + log(0/H) ~ 8.6 (derived from following the ridgeline of 
the data). Thus, empirically calibrated strong-line methods, such as 03N2 and N2, would predict 
oxygen abundances close to this value. However, the galaxies in this range of log(N/0) also extend 
to abundances as low as 12 + log(0/H) ~ 7.9. From a plot of log(N/0) versus stellar mass for 
the SDSS compilation given by Amorfn et al. ( 2010| ) (their Figure 3), our values of log(N/0) are 
roughly 0.6 dex higher than is typical for objects in the same stellar mass range. 



5.5. Expected Oxygen Abundances 

All strong-line calibrations correspond to an assumed electron temperature, so it is informative 
to explore the affects of varying this parameter. Since our MMT spectra include emission lines from 
both 0+ and 0++, we can calculate an expected range in 0/H values based on an assumed reason- 
able possible range in electron temperature. As an example, oxygen abundances were calculated 
for SDSS-082633 assuming electron temperatures of Te = 8 x 10^ K, 1 x 10^ K, and 1.5 x 10"^ K. 
The resulting values of oxygen abundance cover the range of 12 + log(0/H) = 7.5 — 9, and are 
indicated by 'x' symbols in Figure [7} similar results are found for the other 3 galaxies. Thus, we 
cannot rule out any of the various strong line calibration estimates based on this consideration. 

Finally, we can estimate the expected oxygen abundances from the luminosity-metallicity rela- 
tionship and compare these to the strong line estimates in Figure [7j Each calibration is represented 
by a different symbol, and each of our four galaxies is designated with a different color. Using the 
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M-Z relationship of van Zee et al. (2006b 



12 + log(0/H) = 5.67 - 0.151 M^, 



(4) 



we calculate expected oxygen abundances of 8.22, 8.17, 8.07, and 7.97 for the four objects in our 



sample. The M-Z relationship found by Lee et al. (2006b), based on oxygen abundances obtained 



via the direct method and corroborated by Marble et al. ( |2010 ), covers the relevant range in 
luminosity for our sample and is also plotted in Figure [Tj (Note that the range of the T04 relation 
does not extend down to encompass our low- luminosity sample). For each galaxy, a large spread 
is seen amongst the various indicators, highlighting the large uncertainties inherent in metallicity 
determinations . 



6. DISCUSSION 

6.1. Oxygen Abundances from Direct Comparison with Photoionization Models 



The MMT spectra show that three of the objects observed have discrepantly low excitation 
(i.e., low A5007/A3727 ratio; log(032) < —0.4). This implies that the observed nebulae may have 
relatively low ionization parameters and/or be excited by an aging stellar population, which explains 
why we do not see other ions with high ionization potentials (such as He II A4686) in our spectra. 
Since most strong line methods are calibrated for zero-aged main sequence starbursts, and many 
do not incorporate the effects of excitation, one must exercise care in applying these methods to 



determine oxygen abundances (see discussion in van Zee et al. 2006b). In general, the strong line 
methods tend to over-estimate the oxygen abundance for low excitation nebulae, and this offset 
can be as severe as ~0.6 dex ( van Zee et al.|[2006b ). Note that it may be possible to mitigate some 
error by combining empirical and theoretical determinations using an average as [Moustakas et al 



(2010) suggest. 



Additionally, all four of the objects observed have relatively high values of N/0 compared 



to systems of similar stellar mass. As discussed by Yin et al. (2007), PMC09, and Amorm et al 



(2010), the assumption of a single relationship between N/0 and 0/H in photoionization model 
calibrations of the strong line methods can lead to discrepant results when the N abundance does 
not fit this assumed trend. At higher than average values of N/0, the strong line methods tend to 
over-estimate the oxygen abundance. 

The combination of relatively low excitation and relatively high N/0 values leads to a large 
bias in the results from the strong line calibrations. In Figure p[a) we have plotted a comparison of 



the [O II] and [O III] line strengths for our four objects to photoionization models from Stasihska 



Sz Leitherer (1996). Each set of different colored points represents the evolution of a starburst 



'van Zee et al. (2006b I used both direct and empirical oxygen abundances from the literature for 50 objects in 
determining their M-Z relationship. 
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assuming a Salpeter initial mass function (M^p = 100 M©), a stellar mass between 10^-10^ Mq, 
and a metallicity between 12 + log(0/H) = 7.33 and 8.93. Increasing symbol size denotes a 
progression from 1-10 Myr. These models are particularly valuable because they demonstrate the 
additional scatter that can be introduced as the exciting stars of an H II region age, and both 
the ionization parameter and shape of the ionizing spectrum evolve. Figure [Sj^a) shows that H II 
region models with oxygen abundances in the range of 7.9 < 12 + log(0/H) < 8.3 tend to have 
maximum values of R23 (thus defining the turn-around regime in R23; recall log(R23) has an upper 
limit near 1.0). Our four spectra fall among the points for maximum R23 at lower excitation. 
Figure [sj^a) demonstrates that the four objects are inconsistent with relatively high (or relatively 
low) values of 0/H, and, actually, are consistent with the values of 0/H predicted from either the 
luminosity-metallicity relationship or the M-Z relationship. 



Note that in constructing Figure|8j^a), some of the Stasihska &: Leitherer (1996) high-met allicity 
models showed more scatter in this diagram. These models included contributions by Wolf-Rayet 
(W-R) stars to the ionizing spectra, resulting in very hard radiation fields. We have not plotted 
these model points in Figure [sj^a) because our non-detections of He II A4686 indicate very little 
contribution from W-R stars at present. 

MPA-JHU data for SDSS galaxies in the mass range of 7 < log(M^/MQ) < 9 are plotted in 
Figure |8]^b) (gray points). Three of the four objects in the present sample, SDSS-024121, SDSS- 
082639, and SDSS-082633, are discrepantly low excitation nebulae relative to the SDSS sample, 
which tend to have higher excitation. However, most of the P08 "main" sample objects have 
ordinary excitation values. Note that seven of the 24 "main" sample objects did not have z > 0.024 
(consequently, they didn't have [O II] A3727 in their spectra), and so could not be plotted here. 

A complementary analysis is seen in Figure |9] where excitation (P) is plotted against log(R23) 
for this same set of low mass SDSS galaxies (also see |Moustakas et al. ( |2010 ) Figure 12 for a similar 



comparison). As discussed in § 5.3.3, three of our four objects lie well below the mean excitation for 
a given R23 value. Figure [9] also shows an upper limit of log(R23) naturally falls near 1.0, where the 
present sample lies within the turn-around region defined by 0.5 < log(R23) ^ 1-0, corresponding 
to 7.9 < 12 -|- log(0/H) < 8.3. A similar conclusion is drawn from the N2 calibration of PMC09 



(§ 5.3.2) which suggests an oxygen abundance range of 8.30 < 12 -|- log(0/H) < 8.44 for the present 
objects based on the abnormal N/0 variations present (~ 0.1 — 0.2 dex less than other strong-line 
N2 calibrations). While we can't constrain the oxygen abundance more precisely for this sample 



(based on the arguments here in § 6.1) than to put them in the range of 7.9 < 12 + log(0/H) < 



8.4, this is sufficient to disqualify them as high-metallicity M-Z outliers. 

Note that the recalibration of the 03N2 and N2 abundance indicators by PMC09 reduced the 
inferred oxygen abundances, but only by about half of what was needed to bring the objects in line 
with their expected abundances that we are proposing here. Even after recalibration, there is still a 
significant amount of scatter in these relationships. We speculate that taking into consideration the 
additional concern of low excitation would move these points further in the direction of concordance 
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with expectations. 



6.2. The True Nature of the High Oxygen Abundance OutUers 



Although it appears that the "very low mass" , high oxygen abundance galaxies from P08 may 
not be true outliers from the M-Z relation, understanding how these galaxies produced spectra that 
lie far from the median relationships defined by the SDSS galaxies in various diagnostic diagrams 
remains an interesting question. Our MMT observations reveal that these spectra are characterized 
by (1) low excitation (low A5007/A3727 ratios; log(032) < —0.4, see Figures [8] and [9]), (2) relatively 
high N/0 abundance ratios for their stellar masses (see Figure [6]), and (3) relatively deep underlying 
stellar Balmer absorption (ranging from 4 to 8 A). All three characteristics can be explained by a 
previous, but recent burst of star formation. Essentially, these properties are all expected ~ 5 Myr 
or longer after a starburst characterized as a "Wolf-Rayet galaxy." 



Kunth & Sargent ( 1981 ) demonstrated that the presence of strong W-R features in the spec- 



trum of a star forming region are indicative of an intense episode of star formation lasting a relatively 
short (on the order of 10^ years) duration. For an instantaneous starburst, these W-R stars appear 
and disappear during the interval of 3 to 6 Myr ( Leitherer &: Heckman|[T995 ) . Afterward, the most 
massive O stars have all evolved, and the less massive O and B stars produce a significantly softer 



spectrum resulting in lower excitation nebulae. The starburst models of Gonzalez Delgado et al. 



( 1999 ) show that the equivalent width of underlying stellar H/3 absorption increases from roughly 



2.7 to 4.5 A as an instantaneous starburst of solar metallicity ages from to 10 Myr, and continues 
on to nearly 9 A at an age of 100 Myr. On this basis, the MPA-JHU estimates of the Balmer 
absorption are consistent with H II region ages greater than 4 Myr (3.4 A), while our estimates 
(larger due to aperture differences), suggest a lower limit of 7 Myr (5.8 A). Finally, during the 
W-R stage, it is possible for stars to release a large amount of N back into the local ISM (see, e.g., 
[Esteban fc: Vilchez| 1992 ; Kobulnicky et al. 1997 Lopez-Sanchez & Esteban 2010 , and references 
therein), producing inflated N/0 ratios such as is seen in the four MMT spectra. 



6.3. P08 "Main" Sample 

Since all four of our objects, which were picked at random, have roughly identical spectral 
characteristics, we assume that the majority of the "very low mass" sample of P08 are probably 
similar. The "main" sample of P08 was defined by different restrictions, including a redshift cut to 
ensure the coverage of both [O II] and [O III] in the sample spectra (for 17 of the 24 "main" sample 
objects). We obtained the SDSS emission line strengths for the "main" sample from emission line 
analysis by the MPA-JHU collaboration. From these data we found the "main" sample galaxies 



(those with measured [O II] A3727) have normal A3727/A5007 ratios for their masses. In Figure 10 



we show some the properties of the "main" sample relative to the SDSS parent sample, highlighting 
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(a) Oxygen Abundance Models 



(b) Low-Mass SDSS Oxygen Abundances 



Fig. 8. — (a) A comparison of the results from our MMT spectra with results from a large range of 



photoionization models evolved over 10 Myr as determined by Stasinska &: Leitherer (1996) (symbol 
size corresponds to aging sequence). The overlaid dotted contours represent lines of constant R23; 
from left to right log(R23) = 0.0, 0.5, and 1.0. The turn-around region encompasses log(R23) 
values of approximately 0.5 to 1.0. Dashed lines display excitation, where constant excitation of 
[O III]/[0 II] = 1 coincides with the central line, [O III]/[0 II] =0.1 is the low excitation line, and 
[O IIl]/[0 Il]=10 is the high excitation line. Our four spectra are all found near the maximum 
values of R23 for relatively low excitation nebulae, indicating probable 12 + log(0/H) values in the 
range of 7.9 to 8.3. (b) Same plot of [O III] versus [O II] as Figure (a). The numerous gray dots 
are the MPA-JHU data for SDSS galaxies in the range of 7.0 < log(M^/M0) < 9.0, which have 



higher excitation than SDSS-024121, SDSS-082639, and SDSS-082633. Values presented by Yin et 



al. (2007) for a sample of SDSS galaxies containing direct method abundances are also included 



for comparison. Notice that three of our objects from the "very low mass" sample are discrepantly 



low excitation relative to the SDSS and Yin et al. (2007) samples, whereas the "main" sample has 



ordinary excitation values for the most part. 
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Fig. 9. — Excitation parameter P versus the metallicity sensitive R23 index for the four objects 
in this paper (red dots) and the low-mass star- forming galaxies in the SDSS (gray dots). Arrows 
indicate direction of increasing excitation and increasing metallicity for the R23 upper branch. For 
three of these objects, the excitation is rather low for their R23 values. 
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Fig. 10. — Properties of SDSS star-forming galaxies (gray dots) plotted against stellar mass (9 < 
log(M^/M0) < 10). In comparison, the P08 "main" sample (the 17 objects with measurable [O II]) 
are plotted in red. The top panel shows that in contrast to the "very low mass" sample, the "main" 
sample has normal ionization levels. Extremely high N values for the "main" sample relative to 
the SDSS average are shown in the second panel. The bottom two panels show that these objects 
have H/3 absorption equivalent widths and SFRs congruent with the SDSS sample. As all objects 
in the P08 samples exhibit elevated N/0 ratios, we suggest N enrichment is the prominent factor 
influencing these objects. 
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the normal ionization strengths in the top panel, followed by average H/3 absorption equivalent 
widths and SFRs in the bottom two panels. The "main" and 'very low mass" samples are different 
in that the "main" sample galaxies appear to have normal excitations. However, more importantly. 



in Figure 10 we have also plotted N/0 ratio versus stellar mass for the P08 "main" sample, showing 
these galaxies all have very high N/0 values. Similar to the four galaxies in our sample, the P08 
"main" sample data suggest that they have also experienced recent N enrichment. Since high N/0 
values tend to bias strong line methods towards high abundances, and the metallicity offsets for 
the "main" sample are smaller than for the "very low mass" sample, this bias could easily account 
for most of this deviation. 

All objects of the P08 parent sample display the characteristic nitrogen enhancement expected 
for a galaxy after passing through the "Wolf-Rayet galaxy" phase. P08 suggest that these objects 
are "transitional objects running out of fuel": i.e., the high metallicity results because there is only 
a small amount of ISM to enrich. Our overall conclusion is that the high-metallicity outliers are 
special, not due to unusually high oxygen abundances for their mass or rate of star formation, but, 
rather, due to an enrichment of nitrogen relative to oxygen. This nitrogen enrichment may be due 
to being observed during a rare and short lived evolutionary phase. 



CONCLUSION 



We have investigated four of the 41 luminosity-met allicity outliers from the low-mass, high 



oxygen galaxy sample reported by Peeples et al. (2008). From new spectral observations taken at 



the MMT that include the [O II] A3727 line, three of four galaxies were found to have relatively 
low excitation (low A5007/A3727 ratios; log(032) < —0.4) and, more importantly, all four exhibited 
high N/0 values. Each of these characteristics can lead to overestimates of the oxygen abundance 
if standard strong line calibrations are used; furthermore, these two characteristics combined can 
produce the significant metallicity deviations seen for the P08 sample. By comparing our four 



spectra to the photoionization models of Stasihska & Leitherer ( 1996 ) (which include the effects of 
aging of the exciting stars), and the empirical calibrations of PMC09 (which correct for variations 



in the N/0 ratio), we found that significantly lower oxygen abundances are favored. The Stasihska 
fc Leitherer] ( |1996D oxygen abundances are in line with those expected from the luminosities and 
stellar masses of the galaxies, while the PMC09 estimates shift oxygen abundances in the right 
direction. While the "very low mass" sample displays low star formation rates, the normal star 
formation within the "main" sample rules out exhaustive star formation as the variable responsible 
for the unusual sample spectra. With these conclusions, we propose that these galaxies are best 
described as nitrogen enriched, which may result from having recently passed through the "Wolf- 
Rayet galaxy" phase. 
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Table 4. Estimates of the N/0 Ratios for Our Sample 







log(N/0) 








MMT Data 




Temperature 


SDSS-022628 


SDSS-024121 


SDSS-082639 


SDSS-082633 


1.00 X 10'' K 


-0.94 


-1.08 


-0.98 


-1.12 


1.25 X 10" K 


-0.80 


-0.94 


-0.84 


-0.98 


1.50 X 10* K 


-0.70 


-0.84 


-0.74 


-0.88 


Adopted 


-0.80 ± 0.12 


-0.94 ± 0.12 


-0.84 ± 0.12 


-0.98 ± 0.12 



Note. — N/O values for our four galaxy sample, calculated for a range of reasonable 
electron temperatures (12,500 ± 2,500 K). 



